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METHOD AND SYSTEM FOR REDUCING UNDESIRABLE CROSS 
TALK IN DIAGNOSTIC ULTRASOUND ARRAYS 

BACKGROUND 

[0001] Diagnostic Ultrasound is one of the most versatile, lowest cost and 
widely used diagnostic imaging modalities in use today. These systems have a 
controller that controls transmit and receive functions. The conversion of 
electrical energy to acoustic energy and vice versa is performed by a transducer 
array in which individual adjacent elements are separately phased in transmit and 
receive in order to focus the acoustic beam at discrete spatial locations. In the 
ideal ultrasonic array, the adjacent elements operate entirely independently of each 
other. If they do not operate independently - i.e. there is *cross talk' - then the 
beam is distorted and a degraded image results. The cross talk can have its origin 
in either electrical (capacitive) cross talk or acoustic propagation cross talk. 
Generally, electrical cross talk is instantaneous and acoustic cross talk results in a 
delayed cross talk signal due to the much slower propagation velocity of acoustic 
waves versus electromagnetic/electrostatic coupling. Current designs achieve 
cross talk of the order of —30 dB for each of electrical cross talk and acoustic cross 
talk. Further reductions in cross talk may result in image quality improvement. 
Additionally, if an effective means was found for suppressing cross talk then some 
of the intricacies currently employed in transducer design may be avoidable. A 
cost reduction in array design may result. An alternative approach to reducing 
cross talk may also help to facilitate the acceptance of new transducer technologies 
such as electrostatic MEMS transducers that while being very low cost, experience 
significant element to element acoustical cross talk due to their continuous silicon 
substrate. 

BRIEF SUMMARY 

[0002] The present invention is defined by the following claims, and nothing in 
this section should be taken as a limitation on those claims. By way of 
introduction, the preferred embodiments described below include methods and 
systems for reducing undesirable cross-talk in ultrasound arrays. 
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Transmit Mode 

[0003] Each array element is excited by an electrical waveform. Historically 
these waveforms were simple in form ~ impulses, steps or square wave bursts. 
However, the current state of the art is to use a programmable waveform. The 
major driving factor for this is desire to reduce transmitted harmonic signals since 
it is now current practice to form the received signal image at the second harmonic 
of the transmitted signal since the second harmonic signal gives rise to a far higher 
quality image. This imaging technique is typically referred to as "Harmonic 
Imaging". The significance of the need to reduce the transmitted second harmonic 
is that it is preferable to program the shape of the transmitted waveform using 
some form of digital waveform encoding that is fed into a D->A converter prior to 
amplification. Patents describing this technique in detail include Cole U.S. Pat. 
5,675,554 and Hossack U.S. Pat. 5,740,128 (of which are hereby incorporated by 
reference herein). 

[0004] As a central array element is excited by its associated waveform 
generating channel, a smaller signal is generated from its adjacent elements. The 
cross talk signal may be in the range (for example) -20 dB to ^OdB in amplitude 
with respect to the central array element and is delayed by some fraction of a 
microsecond. The actual output of the adjacent element can be measured using a 
needle hydrophone or estimated using a complex transducer model for the 
transducer design in use. PZFlex Finite Element Analysis (Weidlinger Associates, 
Los Altos, CA) is a well known method for comprehensive transducer modeling 
suitable for predicting the cross talk signal. 

[0005] The function of the array can be described by a transfer function matrix 
relating input electrical excitations of one or more elements to measured output 
pressure responses at multiple measurement points in the field. Generally, these 
measurement points are very close to the transducer so that the crosstalk 
component can be more easily distinguished from diffracted components fi-om 
adjacent elements, (i.e. in the very near field, we are measuring the output of the 
adjacent element in isolation from diffracted components from neighbors). 
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[0006] Generally, the initial case with a single element being excited is 
considered with the desire that the output of neighboring elements be as close to 
zero as possible. This condition is satisfied when there is no array crosstalk or 
crosstalk has been adequately suppressed or cancelled. 
[0007] Beamforming generally starts by assuming that each element is 
operating independently and that the total field energy is derived from a linear 
summation of contributions from all elements taking account of diffraction and 
attenuation effects. Thus, the solution of crosstalk problem for the case of a single 
element can be superimposed on the required excitations (and their associated 
required crosstalk cancellation signals) so as to achieve the required beamforming 
effect. Thus, the beamforming task is separated from the crosstalk cancellation 
task. The beamforming may require nonlinear analysis when the sound intensities 
are finite. This is commonly the case in harmonic imaging (imaging at the second 
harmonic of the transmit signal from either tissue or contrast agent bubbles). Even 
in this case, the beamforming process at the array elements is considered in a 
linear fashion - i.e. that all element contributions can be considered in isolation 
and then added when the nonlinear calculation is required. (Typically, the 
nonlinear component is practically zero immediately adjacent to the transducer 
since the harmonic signal accumulates as the wave passes through the media in 
front of the array.) 

[0008] The transfer function is most easily analyzed as an impulse response 
function in the frequency domain. Therefore, in measurements and simulations, 
the time domain responses, H{t), are measured, but are converted to frequency 
domain counterparts, H(J), using a Fourier transform. Conversely, once a required 
excitation is defined in the frequency domain, the time domain equivalent is 
calculated via an inverse Fourier transform. 

[0009] The pressure output. Pi, from the i*'' element can be represented using 
the following matrix equation (Equation 1) where H (abbreviated from H{f)) is the 
transfer function matrix relating input voltages, Vi, to the output pressures, P,. 
This example is for a 9 element array but can be expanded as required. Figure 14 
illustrates this configuration. We will consider the ideal case where an impulse is 



applied to the center element. V5 is an impulse (all 1 's in the frequency domain). 
VI to V4 and V6 to V9 are zero. P5 is the direct output from the 5th element. 
However, since some of the non-diagonal terms in the i/ matrix are non zero, there 
will be finite output from the other elements - i.e. some of Pj to and P<j to Pg 
will be non zero. Once the matrix H is defined, we can solve the matrix problem 
to find values for Vj to such that we do actually get the desired output - i.e. Pj 
to P4 = 0, = ideal impulse response of the 5th element and P^ to Pp =0. 
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[0010] In fact, cross coupling is only significant on a few neighboring 
elements. In the case below, cross coupling is limited to two elements on each 
side of the element being excited. (However, this is not a limitation here — the full 
matrix encompassing all crosstalk terms is considered if desired.) The matrix 
becomes a banded matrix with band width 5. All other values of H are 0. This 
simplifies the characterization of the //matrix and allows for more efficient matrix 
solution. 

[0011] Notice that there is significant symmetry in the matrix. We typically 
assume that crosscoupling (crosstalk) from array element 3 to array element 4 is 
the same as from array element 4 to array element 5 - and so on. There might be 
some irregular behavior at the ends of the array due to different boundary 
conditions. However, even in this case, the change in crosstalk is probably modest 
since it is common practice to have one or two "dummy" elements that are used 
during manufacture to ensure alignment of the saw blade with the required kerf 
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before dicing the 'real' elements. For the matrix shown below, the crosstalk 
function of the two adjacent elements resuhing from impulse excitation of a 
central element is measured. Thereafter, using the 'direct' output from the central 
element and the crosstalk measurements, the banded matrix is populated. 
[0012] For example, we set V5 = 1 (impulse excitation) and Vj to V4 and Vf^ to 
Vg = 0. Hss is defined by the 'direct' output P5. Hs4 and are defined by the 
output on the neighboring element. H^s and H^y are defined by the output on the 
second neighboring element. Thereafter, the symmetry of the matrix is used to 
completely populate the banded matrix. 
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[0013] The impulse responses can be measured experimentally using a small 
hydrophone placed close to each element in a water bath. Suitable hydrophones 
are needle PVDF hydrophones or the 'Golden Lipstick' from Onda Corp, 
Sunnyvale, CA. Alternatively, a prediction of the direct transfer function and 
crosstalk can be made using a suitable computer model — such as using PZFlex 
Finite Element Software (WAX, Los Altos, CA) or ANSYS (ANSYS Inc., 
Canonsburg, PA). Frequently, the true impulse response is not calculated due to 
the practical difficulty of modeling a true impulse. A more practical approach is 
to use a square wave excitation (or similar finite short pulse) of sufficiently high 
bandwidth to produce a relatively unfiltered response. The output response (in the 
frequency domain) can be divided by the frequency domain representation of the 
excitation pulse to derive the impulse response. 
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[0014] The matrix solution can be performed in any of a number of ways. One 
approach is Gaussian ehmination. Another approach is to use the matrix solution 
command in Matlab (Math Works Inc., Natick, MA) (the '\' command, or 
'mldivide'). 

[0015] The approach can be extended to an over constrained case. For 
example, field sample points may be placed closer together in the field than that 
corresponding to array element spacing. For example, one field sample point may 
be placed between array element centers. In this case we consider the example 
where there are 5 array elements and 9 field measurement points. Figure 15 
illustrates this configuration. There are 5 field points lying adjacent to the array 
element centers and 4 points lying intermediate between these 5 field points. In 
this case we have more defined 'output' conditions than defined 'inputs'. Thus, 
the problem is over constrained and needs a pseudoinverse solution. This will find 
a best fitting solution in a least squares sense. The Matlab 'pinv' command 
implements this operation in an efficient and easy to use manner. 
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[0016] In this case, the matrix may also be banded due to finite crosstalk. In 
the following example, crosstalk is limited to the first neighboring element. 
Again, the values of H can be derived either experimentally or using a computer 
model. As before, once a small number of crosstalk measurements are made, the 
matrix can be fully populated by using the symmetry of the system. 
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[0017] In the above two matrix problems, a least squares best fitting solution 
for finding the required excitations of Vj to V3 to produce finite output on only one 
selected element is found using a pseudoinverse operation on matrix H. 
[0018] The above technique can be used in the receive mode. This receive 
mode correction involves (to first order) subtraction of scaled and delayed 
representations of adjacent element signals from the first element's signal. The 
processing to reduce crosstalk is described below in the section: Alternative, non- 
matrix approach to crosstalk reduction calculation. 

[0019] Notice that if a parasitic crosstalk wave is eliminated in either transmit 
or receive, then generally that is sufficient. If the sidelobe resulting from crosstalk 
is eliminated in either transmit or receive, then that sidelobe will be * invisible' in a 
pulse echo operation. 

[0020] The method can be used in all types of medical ultrasound arrays: 
ID, 1.5D, 2D, curved and flat, piezoelectric (PZT, PZN etc,), electrostatic (e.g. 
MEMS or cMUTs) and thermo- acoustic (using a laser to cause transient heating 
over a very small location and acoustic energy generated from the localized 
transient thermal related expansion). 

[0021] Alternative, non-matrix approach to crosstalk reduction calculation: 
[0022] The transfer function of the central transducer is measured by dividing 
(in the frequency domain) the pressure output generated by the associated 
electrical signal input. The observed pressure output from the adjacent element is 
also recorded. Thereafter, it is possible to determine a required electrical signal to 
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be superimposed on the excitation of the adjacent element so as compensate for 
the cross talk in the following manner. The observed cross talk signal (measured 
or simulated) - in the frequency domain - is divided by the transducer transfer 
function as determined above. This calculation determines the electrical excitation 
that is required to give rise to the cross talk signal 'directly' - i.e. not via cross 
talk. If this calculated electrical signal (now in the time domain) is inverted and 
applied to the adjacent element, the direct transduced signal will essentially cancel 
the induced cross talk signal and substantially no net cross talk signal will emanate 
from the adjacent transducer element. The technique can be extended to the next 
adjacent element and so forth. 

[0023] The technique has a recursive element to it in that the electrical signals 
applied to cancel the cross talk will themselves generate a small undesired cross 
talk signal. This second order effect will be largely insignificant where the cross 
talk is modest (e.g. less than -20 dB). If necessary, a recursive waveform design 
approach can be used to properly accommodate these second order cross talk 
effects. As an example, a Monte Carlo approach can be used wherein the shape of 
the excitation waveform is changed in a small way between runs until the best 
overall net cross talk effect is observed taking account of all second, and higher, 
order effects. 

[0024] In an analogous manner, received acoustic waves on a central array 
element will give rise to small, delayed, substantially replica, signals on adjacent 
elements. Therefore, if the cross talk transfer function is known, the anticipated 
cross talk signal on the adjacent element can be subtracted to reveal, at least to a 
good approximation, the original desired cross talk free signal on the adjacent 
element. 

[0025] The following products and services may benefit from the present 
invention: Medical Ultrasound Diagnostics; Non Destructive Examination - 
NDE, NDT; and High Frequency Sonar — Military and civil applications. 
[0026] Exemplary advantages include the following, but are not limited 
thereto: Improved image quality - spatial resolution and contrast quality; 
Potentially, reduced cost; and potentially, improved reliability. 
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[0027] The following are exemplary embodiments of the present invention. 
[0028] 1 . A method and/or system for reducing the impact of undesired inter 
element cross talk in the transmit mode in a diagnostic ultrasound array 
comprising the steps of: 

[0029] a) calculating or measuring the cross talk effect from a first element to 
an adjacent second element; 

[0030] b) applying a first chosen signal to a first element; and 
[0031] c) applying a second chosen signal to a second element wherein the 
excitation to the second element is modified to take account of the calculated or 
measured cross talk effect. 

[0032] 2. A method and/or system for reducing the impact of undesired inter 
element cross talk in the transmit mode in a diagnostic ultrasound array 
comprising the steps of: 

[0033] a) calculating or measuring the cross talk effect from a first element to 
an adjacent second element; 

[0034] b) applying a first chosen signal to a first element; 

[0035] c) applying a second chosen signal to a second element wherein the 

second signal comprises a component that has been added which is a filtered and 

delayed version of the first signal applied to the first element; and 

[0036] d) wherein the filtering and delaying function is a fiinction of the 

calculated, or measured, cross talk effect determined in step a). 

[0037] 3. A method and/or system for reducing the impact of undesired inter 

element cross talk in the receive mode in a diagnostic ultrasound array comprising 

the steps of: 

[0038] a) calculating or measuring the cross talk effect from a first element to 
a second adjacent element; 

[0039] b) receiving a first signal from a first element; 

[0040] c) receiving a second signal from a second element; 

[0041] d) modifying the response of the second element by applying a filter to 

the signal received by the first element and adding the filtered response of the first 

element to the response of the second element; and 
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[0042] e) wherein the fiher is designed to reverse the cross talk effect 
measured in step a). 

BRIEF SUMMARY OF THE DRAWmGS 

[0043] Fig. 1 is an illustration of ideal single element response and degraded 
element response due to cross talk in one embodiment. 

[0044] Fig. 2 is an outline of central array element, with the desired acoustic 
output propagating into tissue in one embodiment. Also shown are the delayed 
and reduced amplitude signals associated with acoustic cross talk and emanating 
from the adjacent elements. 

[0045] Fig. 3 is a block diagram of one embodiment of a process for 

determining electrical excitation of adjacent elements. 

[0046] Fig. 4 is an embodiment of an input spectrum to center element. 

[0047] Fig. 5 shows pulse outputs from central (larger) and adjacent elements 

in one embodiment. The output from the adjacent element is smaller and delayed. 

[0048] Fig. 6 is a spectrum of output of central element in one embodiment. 

[0049] Fig. 7 is an example calculated Transfer Function. 

[0050] Fig. 8 is a spectrum of a cross talk cancellation component. 

[0051] Fig. 9 shows one possible excitation (as a voltage pulse) to an adjacent 

element to compensate for cross talk (dashed line is original excitation to central 

element). As expected, the excitation to the adjacent element is delayed and 

smaller compared to the central element. 

[0052] Fig 10 is a PZFlex™ Finite Element Simulation ~ no compensation 
applied to adjacent element showing Input Gaussian Pulse; Input Pulse Spectrum; 
Output from center element - pulse; Output from center element - spectrum; 
Output from adjacent element - pulse; and Output from adjacent element - 
spectrum. 

[0053] Fig 1 1 is a PZFlex™ Finite Element Simulations - after compensation 
applied to adjacent element showing Input Gaussian Pulse; Input Pulse Spectrum; 
Output from center element - pulse; Output from center element - spectrum; 
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Output from adjacent element — pulse; and Output from adjacent element - 
spectram - Note much smaller signal than for no compensation case. 
[0054] Fig. 12 is one embodiment of a configuration for simulating the cross 
talk only component from a central element to an adjacent element in receive. 
Note the high loss, matched, layer dividing the propagation medium so that the 
adjacent element does not receive a 'direct' signal and that the only signal detected 
is the cross talk signal. 

[0055] Fig. 1 3 is a block diagram of one embodiment indicating the process by 
which the filter to be applied to the adjacent element signal is derived so as to 
remove the contribution due to undesired cross talk. 

[0056] Fig. 14 is an illustration of a 9 element array with 9 field measurement 
points. 

[0057] Fig. 15 is an illustration of a 5 element array with 5 field measurement 
points. 

[0058] Fig. 16 is an illustration of how multiple lasers can be operated in 
phased co-axially aligned. 

[0059] Fig. 1 7 is one embodiment of a modeled transducer. 

[0060] Fig. 18 illustrates one embodiment of a transmit pulse. 

[0061] Fig. 19 shows the pressure on two elements in a response to the pulse of 

Fig. 18. 

[0062] Fig. 20 represents one embodiment of a desired output. 

[0063] Fig. 21 represents cross talk components for application to other 

elements in one embodiment. 

[0064] Fig, 22 represents the outputs responsive to the components of Fig. 21. 

DETAILED DESCRIPTION OF THE INVENTION 
Transmit Mode 

[0065] Each array element is excited by an electrical waveform. Historically 
these waveforms were simple in form - impulses, steps or square wave bursts. 
However, the current state of the art is to use a programmable waveform. The 
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major driving factor for this is the desire to reduce transmitted harmonic signals 
since it is now current practice to form the received signal image at the second 
harmonic of the transmitted signal since the second harmonic signal gives rise to a 
far higher quality image. This imaging technique is typically referred to as Tissue 
Harmonic Imaging'. The significance of the need to reduce the transmitted second 
harmonic is that it is preferable to program the shape of the transmitted waveform 
using some form of digital waveform encoding that is fed into a D->A converter 
prior to amplification. Patents describing this technique in detail include Cole US 
Pat. 5,675,554 and Hossack U.S. Pat 5,740,128. (both of which are hereby 
incorporated by reference herein in their entirety). 

[0066] FIG. 2 shows one system 1 10 for reducing the impact of inter element 
crosstalk in a transmit mode. The system 1 10 is a transmit beamformer 106 and 
associated transducer 112. Additional, different or fewer components may be 
provided, such as including a control processor or path connected with the 
transmit beamformer 106. 

[0067] The transducer 1 12 is an array of piezoelectric or capacitive 
microelectro-mechanical ultrasound transducer (e.g., cMUT) elements 100, 104. 
For microelectromechanical arrays, each element 100, 104 includes a group of 
membranes excited with a same waveform. Different elements 100, 104 of the 
array may be excited with slightly different waveforms, such as a same waveform 
with different timing or waveforms with at least one different frequency or 
amplitude component. The array is a ID, L25D, LSD, 1.75D, 2D, phased, linear, 
curved or other now known or later developed array. 

[0068] The transmit beamformer 1 06 is a transmit beamformer as taught in 
U.S. Patent Nos. 5,675,554 and 5,740,128, but other transmit beamformers with 
different, additional or fewer components may be used. These transmit 
beamformers 106 comprise primarily digital components with digital-to-analog 
converters. Analog transmit beamformers 106 may also be used. 
[0069] The transmit beamformer 1 06 includes a plurality of transmit 
beamformer channels 1 14. The transmit beamformer channels 114 include 
components for generating transmit waveforms including waveform generators 
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108, delays, amplifiers, filters and other now known or later developed 
components. Each channel 114 connects with an element 100, 104 of the array. 
The connection is either permanent or switched, such as with a multiplexer. 
Transmit beamformer channels 1 14 are provided for each element 100, 104 of a 
transmit aperture. Each transmit beamformer channel 1 14 is operable to generate 
a signal or waveform for ultrasound imaging. For example, a one and half cycle 
Gaussian pulse is generated with characteristics adapted for propagation through 
tissue, reflection fi'om tissue boundaries and reception of echoes for imaging. The 
imaging signals are delayed and apodized relative to each other for focusing along 
one or more scan lines. One or more of the transmit beamformer channels 1 14 
also generates a signal for reducing crosstalk from other elements in addition to 
the signal for ultrasound imaging. For example, each transmit beamformer 
channel 114 generates a signal for reducing crosstalk fi'om any adjacent or group 
of adjacent elements 100, 104. 

[0070] Each waveform generator 108 is a memory for generating an entire 
desired waveform with both imaging and crosstalk reduction components, an 
envelope memory, an oscillator, a filter or other now known or later developed 
device for generating a square wave, sinusoidal or other signal at ultrasound 
fi-equencies. Each waveform generator 108 is operable to generate a first 
waveform having a first component to reduce crosstalk fi*om a second waveform 
applied to another element, such as the waveform generated by the waveform 
generator of at least one adjacent element. The transmit beamformer channels 114 
apply the waveforms to the elements substantially simultaneously, such as 
application pursuant to a delay profile for focusing the acoustic energy along a 
scan line in one transmit event. Since the waveforms are applied at a substantially 
same time, the component of one waveform adapted to reduce crosstalk will 
combine with the crosstalk, eliminating or reducing the amplitude of the crosstalk 
firom the adjacent elements. 

[0071] In one embodiment, each waveform generator 108 of each channel 1 14 
generates a single waveform with both a component for reducing the effects of 
crosstalk fi*om adjacent elements and another component for ultrasound imaging. 
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In an alternative embodiment, each channel 114 includes multiple wavefomi 
generators 108 or allows for interleaved generation of waveforms with a same 
waveform generator 108. The crosstalk component is generated separately from 
the imaging component. The separately generated waveforms are then combined 
by one of application to the first element at appropriate timing to provide the 
desired overall waveform or application to a summer in the channel 1 14. 
[0072] hi yet another alternative embodiment, each channel 114 includes one 
or more filters and associated delays connected with adjacent channels 1 14. Each 
filter connects with the channel of another element. The waveform, such as the 
imaging waveform, is filtered by the filter. The passband or filter characteristic of 
the filter is selected as a fimction of the crosstalk response between the elements 
100, 104. The filter may also invert the waveform or otherwise phase modify the 
waveform to make the waveform suitable for crosstalk reduction. Each delay 
delays the filtered waveform, such as to account for differences in focusing delay 
between channels 114 and to align the waveforms for reduction of crosstalk. The 
delayed and filtered waveform from one or more other channels 1 14 is then 
combined with the imaging waveform of the present channel 1 14 as the 
component to cancel or reduce crosstalk. The delay and filtering operation may 
also be combined into a single 'filtering' operation. 

[0073] In a practical imaging system, a group of elements are fired with 
different delays applied to the excitation pulses to each element so as to form a 
focus at a particular point. In this case, the calculated signals to be applied to the 
'adjacent' elements are superimposed (i.e. added ) on to the waveform being used 
on that element to form the focused beam. This waveform addition may be 
performed in a digital domain and stored in the digital waveform memories of a 
multi-channel programmable waveform source - such as shown in US Patent 
5,740,128. Alternatively, the waveforms for particular channels may be added in 
an analog domain for example using well known Operational Amplifier based 
voltage summers. The added or calculated signals reduce the crosstalk by 
canceling out some of the effects of crosstalk. 
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[0074] By removing or reducing crosstalk on transmission, a more desired 
waveform is transmitted. The resulting echoes have less noise or clutter due to 
sidelobes generated by crosstalk. Removal may also be used in the receive 
beamformer, such as by providing filters and delays responsive to crosstalk 
characteristics between elements. Signals form other elements are filtered, 
delayed and added to the received signals to reduce the effects of crosstalk on the 
receive operation. Reduction of crosstalk may be performed on just transmit 
mode, just receive mode or both modes. 

[0075] The received signals are then processed to generate an image with less 
clutter due to crosstalk. This technique can be used in combination with now 
known or later developed diagnostic ultrasound imaging techniques. For example, 
the technique can be used for harmonic imaging modes (all harmonics — fractional, 
sub and super harmonics), with, and without, non-native contrast agents. The 
technique can be used with all forms of Doppler processing — velocity, energy, 
tissue, spectral (Pulsed Wave and Continuous Wave) and all combinations. As yet 
another example, the technique can be used for A-scans, B-scans and C-scans and 
also for 3D. 

[0076] The above system implements a method for reducing the impact of inter 
element cross talk in a transmit mode of a diagnostic ultrasound array. A 
waveform is generated for each element. The waveforms include components 
operable to reduce crosstalk from application of other waveforms at other 
elements. As discussed above, the components for reducing crosstalk are 
generated with imaging waveforms as a single waveform, such as by modifying 
the desired imaging pulse to take account of crosstalk from one or more other 
elements. The components may alternatively be generated by filtering and 
delaying the desired imaging pulses from other channels 1 14 as a function of a 
crosstalk transfer function. In yet another alternative, the desired imaging pulse 
and any crosstalk reduction components are generated separately and combined. 
The two components are combined by applying both to the same element with 
desired relative timing or summing prior to application to the element. 



-16- 



[0077] Figure 4 shows the spectrum of a desired imaging signal or pulse for 
application to one or more elements. Figure 5 shows the output (solid line) from 
an element of a time domain waveform corresponding to the spectrum of Figure 4. 
The dashed line waveform of Figure 5 shows the crosstalk pressure response from 
application of the same waveform to an adjacent element. The crosstalk pressure 
is delayed and has a smaller amplitude. Figure 6 shows the spectral result of the 
desired waveform and the undesired crosstalk. Figure 7 shows a calculated 
transfer function. The transfer function represents the difference in spectral 
response from the desired response of Figure 4 due to the undesired crosstalk (i.e., 
the spectrum of Figure 6 minus the spectrum of Figure 4). Figure 8 shows the a 
spectral response of a component waveform for reducing the effects of the 
crosstalk. Figure 9 shows the imaging pulse or signal (dashed line) with the time 
domain waveform or component for canceling crosstalk from the adjacent 
element. The crosstalk reduction component is delayed and of smaller amplitude 
relative to the desired imaging pulse. The crosstalk reduction component is 
substantially an inverse waveform to the crosstalk generated on the element from 
the adjacent element. 

[0078] As further examples. Figures 10 and 1 1 show simulated crosstalk. 
Figure 10 shows an input Gaussian pulse and the associated spectrum applied to 
an element. The output of the element and the associated spectrum are also 
shown. Due to the oscillation created in the element, the ringdown time of the 
input waveform is extended. The crosstalk generated at an adjacent element and 
the associated spectrum is also shown. Figure 1 1 shows the same information, but 
with the component to reduce crosstalk added with the crosstalk at the adjacent 
element. The pressure scale for the cross talk is much lower. The output 
waveform of the primary element (e.g., center element) is different in Figures 10 
and 1 1 since the crosstalk cancelation signal causes the output of the central 
element to change slightly since the cancellation signal is couple back to the 
primary element. In alternative embodiments, the cancellation or crosstalk 
reduction component is added to the waveform of one element to reduce the 
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crosstalk at a different element, such as adding the component to the Gaussian 
pulse shown in Figures 10 and 11. 

[0079] Desired imaging pulses, such as the Gaussian pulses shown in Figures 
10 and 1 1, are applied to each element in the transmit aperture. Any components 
for reducing crosstalk from adjacent or other elements are also applied. The 
waveforms having both the desired imaging pulse and the crosstalk reduction 
pulse are applied to all the elements substantially simultaneously. Substantially is 
used to account for the relative delays between channels 1 14 for focusing as well 
as delays in the crosstalk component to account for propagation of the crosstalk 
energy between elements. The crosstalk component is operable to reduce the 
crosstalk from all elements, a subset of elements, adjacent elements or groups of 
adjacent elements. 

[0080] A similar process is used for reducing the impact of inter element 
crosstalk in a receive mode in a diagnostic ultrasound array. Echo signals are 
received at each element across a receive aperture of the array. To reduce the 
crosstalk from reception of echo signals at other elements, the received signals are 
modified as a fiinction of the crosstalk characteristic between the elements. For 
example, the received signal at one element is filtered as a fiinction of a crosstalk 
characteristic, such as a crosstalk transfer fiinction between the elements. The 
filtered signal represents an inverse of the likely crosstalk component at another 
element. The filtered signal is then added to the received echo signal at the other 
element, reducing the effects of any crosstalk. The crosstalk transfer fiinction or 
effects of the crosstalk are measured or calculated as discussed below. 
[0081] Referring to Fig. 1, as a central array element 100 is excited by its 
associated waveform generating channel, a smaller signal 102 is generate from its 
adjacent elements 104. The cross talk signal may be in the range (for example) - 
20 dB to — 40dB in amplitude with respect to the central array element and is 
delayed by some fraction of a microsecond. The actual output of the adjacent 
element can be measured using a needle hydrophone or estimated using a complex 
transducer model for the transducer design in use. PZFlex Finite Element 
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Analysis is a well known method for comprehensive transducer modeling suitable 
for predicting the cross talk signal. 

[0082] 'Central' and 'adjacent' are used for convenience. The technique is 
applicable to any or all array elements except the end elements for which there is 
no 'adjacent' element on the outer edge. The 'central' element 100 is the primary 
desired element being driven (first order effect). The adjacent element 104 
vibrates - in addition to the vibration due to its own primary (1st order) voltage 
drive, according to a second order effect that is due to primary effect on its 
neighbor — the central element. Fig 2 shows another representation of the this 
undesired cross talk. 

[0083] The function of the array can be described by a transfer function matrix 
relating input electrical excitations of one or more elements to measured output 
pressure responses at multiple measurement points in the field. Generally, these 
measurement points are very close to the transducer so that the crosstalk 
component can be more easily distinguished from diffracted components from 
adjacent elements, (i.e. in the very near field, the output of the adjacent element is 
measured in isolation from diffracted components from neighbors). 
[0084] Generally, the initial case has a single element being excited as shown 
at act 302 in Fig. 3. The output of neighboring elements is desired to be as close 
to zero as possible. This condition is satisfied when there is no array crosstalk or 
crosstalk has been adequately suppressed or cancelled. After applying a known 
excitation to the central element 100, the output at the central element 100 and 
adjacent elements 104 are measured in act 304. In act 306, a fast Fourier 
transform is performed on the electrical excitation applied to the central element in 
act 302 as well as the output of the central element measured in act 304. In act 
308, the central element transfer function is calculated from the transformed 
information. The transfer function is the ratio of output pressure waveform (in the 
frequency domain) to input electrical signal (in the frequency domain). In act 310, 
the fast Fourier transform of the measured signal at the adjacent element is divided 
by the transfer function calculated in act 308. In act 312, the result of act 310 is 
inverted. The signal is converted to the time domain via inverse Fourier 
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transform, inverted to produce the quantity which when added to the original 
crosstalk sums to zero, and then converted back to the frequency domain. 
However, as a practical matter, in act 3 12, the result of act 3 10 may be inverted 
directly in the frequency domain - i.e. phase rotated by 1 80 degrees at every 
frequency sample. The result of act 3 12 is then divided by the transducer transfer 
function to determine the first order approximation to the required additional 
electrical stimulation to the adjacent element to cause crosstalk cancellation. The 
result of this division is then converted back to the time domain via an inverse 
Fourier transform and is then superimposed on the signal applied to the adjacent 
element to result in reduced crosstalk. 

[0085] The above process demonstrates one example of determining a 
crosstalk characteristic from one element to another second element where the 
elements are adjacent to each other in the array. The same process may be used 
for non-adjacent elements. The crosstalk characteristic in the process of Figure 3 
relies on measuring the crosstalk. Alternatively, the crosstalk is modeled or 
calculated. The crosstalk reduction component is generated as a function of a 
characteristic of the crosstalk, such as the amplitude, delay and frequency of the 
crosstalk. 

[0086] Figure 13 shows a process similar to Figure 3 where a filter response is 
determined as a function of the crosstalk characteristic. The example of Figure 13 
is for a receive mode filter, but the process may be altered for application to a 
transmit mode filter by using transmitted signals. In act 1302, a known acoustic 
excitation is applied to a central element and a responsive echo signal is received. 
Alternatively, a likely receive signal is applied to the central element. In act 1304, 
the output voltage from the central and an adjacent element responsive to act 1302 
is measured. The outputs voltage is the received signal generated in response to 
transmission using the known excitation. Alternatively, the output is calculated or 
simulated. In act 1 306, a fast Fourier transform is performed on the electrical 
received signals on the central element and adjacent elements measured in act 
1304. In act 1308, the central element to adjacent element transfer function is 
calculated from the transformed information. This central element to adjacent 



-20- 



element transfer function is calculated by dividing the electrical output from the 
adjacent element (in frequency domain) by the pressure wave input to the central 
element (in frequency domain). The 'inversion' act may occur in one of several 
places. In act 1310, the transfer fiinction calculated in act 1308 is inverse fast 
Fourier transformed. The polarity of the transfer function is reversed either prior 
to or afler the inverse transform. The result represents the desired filter response 
for reducing crosstalk. In act 1312, a received signal used for imaging is applied 
to a filter programmed with the desired filter response. For example, a received 
signal from the adjacent element is applied to the filter. The filtered output is 
added to the received signal at the center element to reduce any crosstalk. This 
process may be repeated for designing the filters of additional, all or a subset of 
elements of the receive aperture. Multiple filters may be used for any given 
channel to reduce crosstalk from a plurality of elements in any given element. 
[0087] For either the transmit or receive modes, the measurements are 
performed for every element. Alternatively, the symmetry of the array is used. In 
one embodiment, a banded matrix is populated as function of the crosstalk 
characteristic measured at one, two or other subset of elements. The resulting 
banded matrix is used to determine the component waveform for reducing 
crosstalk at each element due to transmit or receive signals at one or more other 
elements. 

[0088] Beamforming generally starts by assuming that each element is 
operating independently and that the total field energy is derived from a linear 
summation of contributions from all elements taking account of diffraction and 
attenuation effects. However, cross-talk occurs so that a cross-talk signal is 
measured in act 304. Thus, by solving the crosstalk problem for the case of a 
single element, then the excitations and their associated crosstalk cancellation 
signals are superimposed so as to achieve the required beamforming effect, 
reduction or elimination of cross talk across the array. The beamforming task is 
separate from the crosstalk cancellation task. The beamforming may require 
nonlinear analysis when the sound intensities are finite. This is commonly the 
case in harmonic imaging (imaging at the second harmonic of the transmit signal 
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from either tissue or contrast agent bubbles). Even in this case, the beamforming 
process at the array elements is considered in a linear fashion — i.e. that all element 
contributions can be considered in isolation and then added when the nonlinear 
calculation is required. Typically, the nonlinear component is practically zero 
immediately adjacent to the transducer since the harmonic signal accumulates as 
the wave passes through the media in front of the array. 

[0089] The transfer function is most easily analyzed as an impulse response 
function in the frequency domain. Therefore, in measurements and simulations, 
the time domain responses, H(t)^ is measured, and then converted to the frequency 
domain counterparts, H(f)y using a Fourier transform. Conversely, once a required 
excitation is defined in the frequency domain, the time domain equivalent is 
calculated via an inverse Fourier transform. 

[0090] The pressure output, P„ from the i^** element can be represented using 
the following matrix equation (Equation 1) where H (abbreviated from H(f)) is the 
transfer function matrix relating input voltages, F,, to the output pressures, P,. 
This example is for a 9 element array such as shovm in Figure 14, but may be 
expanded as required. An impulse is applied to the center element so that is an 
impulse (all 1 's in the frequency domain). V] to V4 and to Fp are zero. Ps is the 
direct output from the 5th element. However, since some of the non diagonal 
terms in the H matrix are non zero, there will be finite output from the other 
elements - i.e. some of Pj to P4 and P^ to Pp may be non zero. Once the matrix H 
is defined, the matrix problem is solved to find values for Vj to Vg such that the 
desired output - i.e. Pj to P4 = 0, P5 = ideal impulse response of the 5th element 
and P^ to Pg =0 results or such that at least some cross talk is reduced for one or 
more elements. 
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[0091] In fact, cross coupling may only be significant on a few neighboring 
elements. In the case below, cross coupling is limited to two elements on each 
side of the element being excited. The matrix becomes a banded matrix with band 
width 5. All other values of H are 0. This simplifies the characterization of the H 
matrix and allows for a more efficient matrix solution. 
[0092] Notice that there is significant symmetry in the matrix. The 
crosscoupling (crosstalk) from array element 3 to array element 4 may be the same 
as from array element 4 to array element 5 — and so on. There might be some 
irregular behavior at the ends of the array due to different boundary conditions. 
However, even in this case, the change in crosstalk may be modest since it is 
common practice to have one or two *dummy' elements that are used during 
manufacture to ensure alignment of the saw blade with the required kerf before 
dicing the *rear elements. For the matrix shown below, the crosstalk function of 
the two adjacent elements resulting from impulse excitation of a central element is 
measured without measuring at other elements. Thereafter, using the 'direct' 
output from the central element and the crosstalk measurements, the banded 
matrix is populated. 

[0093] For example, Fj = 1 (impulse excitation) and Vj to V4 and K<y to Fp = 0. 
Hss is defined by the ^direct' output Pj. and H^^ are defined by the output on 
the neighboring element. H^s and i/57 are defined by the output on the second 
neighboring element. Thereafter, the symmetry of the matrix is used to 
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completely populate the banded matrix, providing cancellation wave forms for 
each element Vj-Vg. 
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[0094] The impulse responses can be measured experimentally using a small 
hydrophone placed close to each element in a water bath. Suitable hydrophones 
are needle PVDF hydrophones or the 'Golden Lipstick' from Onda Corp., 
Sunnyvale, CA. Altematively, a prediction of the direct transfer function and 
crosstalk can be made using a suitable computer model - such as using PZFlex 
Finite Element Software (WAI, Los Altos, CA) or ANSYS (ANSYS Inc., 
Canonsburg, PA). Frequently the true impulse response is not calculated due to 
the practical difficulty of modeling a true impulse. A more practical approach is 
to use a square wave excitation (or similar finite short pulse) of sufficiently high 
bandwidth to produce a relatively unfiltered response. The output response (in the 
frequency domain) can be divided by the frequency domain representation of the 
excitation pulse to derive the impulse response. 

[0095] The matrix solution can be performed in any of a number of ways. One 
approach to the solution is Gaussian elimination. Another approach is to use the 
matrix solution command in Matlab (the '\' command, or 'mldivide'). 
[0096] The approach can be extended to an over constrained case. For 
example, field sample points may be placed closer together in the field than that 
corresponding to array element spacing. For example, one field sample point may 
be placed between array element centers. For example, 5 array elements and 9 
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field measurement points are used. Figure 15 illustrates this configuration. There 
are 5 field points lying adjacent to the array element centers and 4 points lying 
intermediate between these 5 field points. In this case we have more defined 
'output' conditions than defined 'inputs'. Thus, the problem is over constrained 
and needs a pseudoinverse solution. A best fitting solution is found by a least 
squares analysis. The Matlab 'pinv' command implements this operation in an 
efficient and easy to use manner. 
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[0097] In this case, the matrix may also be banded due to finite crosstalk. In 
the following example, crosstalk is limited to the first neighboring element. Again, 
the values of H can be derived either experimentally or using a computer model. 
As before, once a small number of crosstalk measurements are made, the matrix 
can be fiiUy populated by using the symmetry of the system. 
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[0098] In the above two matrix problems, a least squares best fitting solution 
for finding the required excitations of VI to V5 to produce finite output on only 
one selected element is found using a pseudoinverse operation on matrix H. 
[0099] If a parasitic crosstalk wave is eliminated in either transmit or receive, 
then generally that is sufficient. If the sidelobe resulting from crosstalk is 
eliminated in either transmit or receive, then that sidelobe may be 'invisible' or 
have minimal undesired effect in a pulse echo operation. ^ 

Receive Mode 

[00100] The above technique can also be used in the receive mode. The 
correction involves (to first order) subtraction of scaled and delayed 
representations of adjacent element signals from the first element's signal. 
[00101] Received acoustic waves on a central array element vsdll give rise to 
small, delayed, substantially replica, signals on adjacent elements. Therefore, if 
the cross talk transfer function is known, the anticipated cross talk signal on the 
adjacent element can be subtracted to reveal, at least to a good approximation, the 
original desired cross talk free signal on the adjacent element. 
[00102] In the receive mode, the transfer function relating the echo received 
on the adjacent element to the echo received on the central element is calculated. 
In a practical situation^ the 'real' signal on the adjacent element is compounded 
with the unwanted cross talk signal from the central element. Therefore, this 
calculation of transfer function is best performed using a detailed simulation - e.g. 
PZFlex finite element analysis (Weidlinger Associates Inc. Los Altos, CA). In 
this simulation, the artificial constraint is applied that the input pressure function is 
received by the central element but not by the adjacent element. This may be 
realized by dividing the acoustic medium through which the source pressure 
function is propagating into two with an extremely high loss, acoustically 
matched, medium as a dividing layer down the model to prevent any 'real' input 
being received by the adjacent element. Fig. 12 shows this configuration. This 
model produces the desired acoustic output from the central element (time 
domain) and the undesired cross talk signal on the adjacent element (also time 
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domain). The FFT of both of these functions is calculated and the response 
(frequency domain) of the adjacent element divided by the response (frequency 
domain) of the central element. This defines the transfer function of the cross talk 
energy (frequency domain). A digital filter that calculates the signal to be added 
to adjacent element so as to cancel just the cross talk signal may be found by 
taking the IFFT of the transfer function and applying a change of sign - i.e. a 
negative prefix to the digital filter coefficients. Thereafter, in a practical system, 
the cross talk response may be removed by convolving this filter with the signal 
from the central element and adding it to the response of the adjacent element to 
)deld the cross talk free adjacent element signal. This determination is repeated 
for each element in the receive aperture. 

[00103] Of course, for reasons of economy, a more simplistic approach may 
be desired. As one example, a scaled (up or down) and delayed and/or phased 
replica of the center element response is added, or subtracted, from the response of 
the adjacent element to yield an approximate cross talk fi-ee response for the 
adjacent element. 

[00104] This derivation solves only for the first order correction. In 
practice, the second order effects, such as the cross talk effect of the adjacent 
element onto the central element - before the signal to be input to the cross talk 
compensating filter - has not been accounted for. The basic technique can be 
extended to take account of these second (and higher) cross talk effects by more 
complex analysis — including, for example, a Monte Carlo approach to finding the 
required filter coefficients that properly remove the first, second and higher cross 
talk effects. In practice, it is believed that the first order correction descfribed 
herein is sufficient. 

Variations 

[00105] The method is particularly well suited to use with transducers that 
have inherently high levels of cross talk. In particular, silicon MEMS transducers, 
typically exhibit significant cross talk. These MEMS transducers are described in 
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detail in the patents by Khuri-Yakub et al. These patents (hsted below) are hereby 

incorporated by reference herein in their entirety: 

[00106] 1. 5,6 1 9,476 'Electrostatic ultrasonic transducer' 

[00107] 2. 5,870,351 'Broadband microfabriated ultrasonic transducer and 

method of fabrication' 

[00108] 3. 5,894,452 'Microfabricated ultrasonic immersion transducer' 
[00109] 4. 5,982,709 'Acoustic transducers and method of microfabrication' 
[00110] 5. 6,004,832 'Method of fabricating an electrostatic ultrasonic 
transducer' 

[00111] The filtering operations may operate on the raw radio frequency 
signals. Altematively, the processing may operate on complex baseband 
demodulated signals as are common in modem ultrasound systems. The 
processing may also occur in the analog domain using analog delay components, 
analog filters and analog summers. 

[00112] The method also has applicability in the context of laser generated 
ultrasound. See the following references: 

[00113] Buma T, Spisar M, O'Donnell M. 'High-fi"equency ultrasound array 
element using thermoelastic expansion in an 

elastomeric film'. Applied Physics Letters, vol.79, no.4, 23 July 2001, pp.548-50. 
Publisher: AIP, USA 

[00114] Buma T, Spisar M, O'Donnell M. 'Thermoelastic generation of 
ultrasound using an erbium doped fiber amplifier' 1999 IEEE Ultrasonics 
Symposium. Proceedings. International Symposium (Cat. No.99CH37027). 
IEEE. Part vol.2, 1999, pp.1253-6 vol.2. Piscataway, NJ, USA 
[00115] In these cases, a laser is fired at a surface comprising a material with 
a high coefficient of thermal expansion (e.g.. Chromium or PDMS). The rapid 
expansion caused by the laser pulse gives rise to an acoustic wave. A second, or 
more, laser is fired at a point on the surface adjacent to the first point. The 
excitation of this second laser is delayed and attenuated as described in the prior 
discussion so that the resulting acoustic wave due to the second laser at the second 
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point cancels the result of the cross coupled wave emanating across the surface 
from the first point. 

[00116] Multiple co-axially aligned laser spots are used with successive 
lasers focused to wider and wider spots. The center of the spot receives 
contributions from all lasers. Successive ^rings' going outwards receive 
contributions from fewer laser spots. Therefore, the approach taken is to define 
the excitation functions desired for each of the rings. In Figure 16, 3 rings from 3 
lasers are shown. The excitation function required for the outer ring is applied to 
Laser 3. The Laser 2 excitation is calculated by taking the desired function for 
ring 2 and subtracting from this the contribution already being made by Laser 3. 
Finally, the contributions from Laser 2 and Laser 1 are subtracted from the desired 
fiinction for ring 1, and this residual signal is applied to Laser 1 . The excitation 
functions may be approximations since only 'positive' pressures can be generated 
- since only heating and associated expansion is available. The idea here would 
be to superimpose the outer rings so that the waves generated by the heating in the 
outer rings causes an acoustic wave approximately 180 degrees out of phase with 
respect to the wave propagating outwards from the center spot. 
[00117] The Matlab code below is operable to read output from PZFlex for 
input excitation and output from central and adjacent elements and then calculate 
the required additional excitation to be applied to an adjacent element so as to 
cancel the effect of cross coupling. 



[00118] 


% Matlab program to do the decoupling calculation 


[00119] 


% Read PZFlex finite element data files for 


[00120] 


% 1 Input to 'central' element — assume Flex files are in matlab form 


[00121] 


% 2 Output from 'central' element 


[00122] 


% 3 Output from adjacent element 


[00123] 


load input 


[00124] 


load output_center 


[00125] 


load output_adjacent 


[00126] 


load time 


[00127] 


ns=8 1 92; % Use 2''N samples 
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[00128] ts=time(2)-time(l); % Time axis 

[00129] t=(l:ns)*ts; 

[00130] input_fft=fft(input); % FFTs 

[00131] output_center_fft=fR(output_center); 

[00132] output adjacent_fft=fft(output adjacent); 

[00133] fTnax=(l/ts)/2; % Maximum frequency 

[00134] fs=fmax/ns; % Sample frequency 

[00135] f=(l :ns)*fs; % Frequency axis 

[00136] % Transfer fiinction calculation 

[00137] trans_fiinc=output_center_fft./input_fft; 

[00138] % required excitation for the input to the adjacent element 

[00139] reqd excit=output_adjacent_fft./trans_fiinc; 

[00140] reqdj)ulse=-real(ifft(reqd_excit)); % Required excitation to 

adjacent 

A. Transmitting Mode Model 

[00141] A 2D ultrasound transducer array model is modeled using PZFlex. 
The model is illustrated in Figure 17. The actual FEA model takes advantage of 
the symmetry plane and hence only the FEA models one half of the device 
illustrated. Important acoustic properties and geometry design parameters are 
described in the Tables 1 and 2. 



Table 1. Acoustic properties for the FEA model 



Material 


Impedance (MRayl) 


PZT 


34.2 


Matching material 


3.92 


Backing material 


2.15 


Lens 


1.5 



Table 2. Geometry design parameters for the FEA model 



Composite layer thickness 


0.70mm 


Element pitch 


0.30mm 


Backing layer thickness 


1.00mm 


Elevation Length 


10.0mm 


Matching layer thickness 


0.35mm 


Kerf width 


0.10mm 


Lens layer thickness 


0.30mm 


Kerf cut depth into backing 


0.08mm 



-30- 



[00142] There are seven elements indexed from 1 to 7 simulated in this 
model, where 4 is the center element (Figure 4). The pressure output signals are 
measured at seven field points immediately above the center of elements. In these 
simulations we assume a 50 Q source impedance. When we transmit on the 
central element, we set the voltage to zero at the source for the adjacent elements - 
so that the adjacent elements are properly terminated with 50 CI to ground. The 
adjacent elements are not directly shorted to ground. This approach is necessary 
in order to apply the current superposition theorem correctly. The circuit is 
illustrated in Figure 5. We use a wideband Gaussian pulse (2.5MHz, 80% -6dB 
fractional bandwidth) as the central element excitation. Time domain records of 
pressure outputs were transferred to a PC mnning Matlab (Mathworks, Natick, 
MA). Using an FFT in Matlab, we converted the time domain records to their 
equivalent frequency domain responses and form the frequency domain transfer 
functions and crosstalk functions. These responses were then used to populate the 
matrix equation (Equation 1). Our initial calculations use as many pressure 
measuring locations as transducer source points. Therefore, we are using a square 
matrix amenable to standard Gaussian elimination. 

[00143] Figure 18 illustrates the applied input pulse to (center) Element 4 in 
both time and frequency domain while the other elements are grounded. Figure 19 
illustrates the pressure outputs from Elements 4 and 5, respectively, in both time 
and frequency domain. It is evident from these figures that the crosstalk level in 
the transmit mode is approximately -18 dB. To first order, the crosstalk signal is a 
delayed and attenuation version of the desired signal output from Element 4, 
[00144] Figure 20 illustrates the selected desired output waveform from the 
central element. In the interests of transmission efficiency, this waveform is 
chosen to be somewhat similar to the observed output for a wideband Gaussian 
excitation. Figure 21 illustrates the excitation fiinctions to be applied to each of 
elements 1, 2, 3 and 4, where the 4th is the central element. The functions applied 
to 5, 6 and 7 correspond to those applied to 3, 2 and 1 respectively due to the 
symmetry of the structure. These functions were derived from the solution to 
Equation 1 . The actual matrix solution is performed in the frequency domain and 
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the waveforms illustrated in Figure 22 are obtained via an inverse Fourier 
transform. 

[00145] Figure 22 illustrates the pressure output obtained from the central 
and adjacent elements when the excitation functions illustrated in Figure 21 are 
applied to the transducer terminals. A significant reduction in crosstalk has been 
achieved. Crosstalk has been reduced from -18dB to -27 dB. 
[00146] The crosstalk reduction process is repeated using approximately 
twice as many pressure field sample points - one in front of the center of each 
element and one in additional point midway between these sample points - i.e. 
above the intermediate kerf region. In this case there are 13 field points for the 7 
transducer elements. In fact, due to symmetry, only 7 field points and 4 transducer 
elements are analyzed. Again, a broadband Gaussian pulse, identical to that used 
above, is used as the desired central element output. The excitation functions are 
obtained afler finding a pseudo-inverse solution to the matrix equation (Equation 
3). The actual waveforms calculated using this second approach, and the 
consequent reduction in crosstalk, are similar to those for the square matrix case 
(Equation 1). 

[00147] In one embodiment, an array of waveforms are generated for 
forming a focused beam. The waveform applied to each or a subset of elements 
has a component that is a delayed, filtered and possibly inverted form of the 
excitation of an adjacent element. The beamformed imaging signals have 
differential delays to form the focus but the crosstalk canceling signals have 
uniform delays with respect to immediate neighbors since propagation time 
between adjacent elements is always the same. 

[00148] While the invention has been described above by reference to 
various embodiments, it should be understood that many changes and 
modifications can be made without departing from the scope of the invention. It is 
therefore intended that the foregoing detailed description be regarded as 
illustrative rather than limiting, and that it be understood that it is the following 
claims, including all equivalents, that are intended to define the spirit and scope of 
this invention. 



